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ABSTRACT 

We discuss approximation techniques for use in 
numerical schemes for estimating spatially varying 
coefficients In continuum models such as those for 
Euler-Bemoulll beams. The techniques are based on 
qulntlc spline state approximations and cubic spline 
parameter approximations. Both theoretical and 
numerical results are presented. 

Introduction 

We discuss here our efforts on the development of 
numerical algorithms for the estimation of parameters 
in variable structure elastic models. The class of 
problems we Investigate Is of fundamental Importance 
in the use of continuum models (such as those for 
beams, plates, thin shells, etc.) to represent large 
flexible space structures [1-4]. Our own research has 
been motivated by the need to detect and/or estimate 
structural / material property changes in such struc- 
tures while they are on orbit. In this regard the 
usefulness of parameter estimation techniques (l.e., 
inverse algorithms) for spatially varying parameters 
in models such as 

p(x) (EI(x) ) * f(6,t,x) (1) 

3t c 3x‘ tx 

from the Euler-Bemoulli theory or In more sophisti- 
cated models arising in the Timoshenko theories should 
be rather obvious. In this presentation we describe 
spline-based methods that use and extend In a non- 
trivial way the Ideas developed in [5-7] for constant 
parameter systems. These extensions allow one to treat 
variable parameters (such as mass density p, elastic 
modulus El, load parameters s) In models such as (1) 
and were first developed in the context of variable 
coefficient first order parabolic systems in [8]. 

Although our Ideas have much wider applicability 
in elastic structures, to illustrate the basic princi- 
ples we, for simplicity, restrict our considerations 
here to a simply supported beam with constant mass 
density and variable stiffness, which Is under a known 
load f. We assume that we have normalized our model 

and thus consider (with D * 

u tt ♦ 0 2 (a(x)D 2 u) * f(t,x) 0<x< 1, t»0, (2) 

with boundary conditions 

u(t,0) * u ( t , 1 ) * u xx (t,0) * u xx (t,l) * 0 (3) 

and initial conditions 

u(0,x) - 4 ( x) 

(4) 

u t (0,x) * *(x) . 

The methods presented here readily extend to more 
general cases (e.g., unknown parameters in more general 
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boundary conditions. In Initial conditions or load, 
variable mass and damping, etc.) as will be discussed 
In a forthcoming paper. 

The problem we consider is that of estimating 
functions o In (2) from observations of the state u. 

To be more specific, we consider a least squares flt- 
to-data for (2). That Is, given observations y^j for 

for uft^Xj;®), where u Is the solution to (2)-(4) 

corresponding to the parameter a from the set Q of 
admissible parameters, we seek a parameter a* that 
minimizes 

j(a) - I |u(t 1 f Xj;a) - y^] 2 (5) 

1 , j " 

over Q. 

We outline briefly the Ideas of our approach to 
this problem. We first rewrite (2)-(4) as an abstract 
system In an appropriately chosen Hilbert space Z * 
Z(a), our state space. The resulting abstract Iden- 
tification problem Is then approximated on finite 

u 

dimensional state subspaces Z (a) which In our case 
here are generated by qulntlc spline elements. We 

M 

then obtain estimates a from minimizing a fit-to-data 
criterion using the approximate states. 

u 

The problem of obtaining a Is, however, not 
computationally feasible since It requires that a min- 
imization procedure be carried out over an Infinite 
dimensional function space or set Q. We therefore M 
Introduce a further approximation by emptying sets Q n 
(In the results reported below these Involved cubic 
spline Interpolations to the elements In Q) In the 
minimization problems. This results In a double 
approximation procedure: a state approximation (spaces 

Z to approximate Z) and a parameter approximation 
M 

(sets Q to approximate Q); thus computationally we 

N M 

obtain least squares parameters a ’ and any conver- 
gence discussions must Involve a double limit process 
for N M + •. As we shall note below, the Trotter- 
Kato approximation theorem for linear semigroups, which 
we have successfully used In a t.jmber of other problems 
[5, 6, 7 ,9], can be used In the present case to establish 
a convergence theory. (For the convenience of the 
reader, we note that we are using the notation a for 
the unknown parameter function here whereas In the 
general formulations of [5-9] the symbol q Is used to 
denote a vector of unknown parameters.) We proceed to 
outline the schemes we have used, discuss briefly con- 
vergence results, and present some of our numerical 
findings. 

Qulntlc/Cubic Spline Approximation Schemes 

Given (2)- (4), we rewrite this system In abstract 
form In the parameter dependent state space 

Z(o) 5 (H 2 0 Hq) x H^. Here H 2 Is the space of 
elements in H 2 (0,1) with Inner product 
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< 2 ,w> a ■ <a0 2 *,,0 2 w 1 > 0 * <i 2 .*» 2 > 0 where < , > Q Is the 

usual Inner product In H®(0,1) ■ LjlO.l). Define the 
operator A(a) by 

/ 0 1 

A(a) » [ 

\ -D 2 (aD 2 ) 0 

on dom(A) ■ (i£ H 4 * H 4 |zj(n) ■ z 2 f 1 ') “ ^ 2 z^(n) * 
0 2 z 2 (r.) ■ 0. n ■ 0,1). Then putting z« (u,u t ) T , 

Zq * ( * . * ) T . and F » (f,0) T , we can rewrite the system 
( 2)- ( 4) as 

z(t) * A(a)z(t) ♦ F(t) t > 0, 

( 6 ) 

Z( 0 ) * Zg. 


The estimation problem is then one of minimizing 


J(») - Ii 1 (t 1 ;«)(x j ) - ‘*.j| 2 (7) 

2 

over Q, where we assume Q Is a given subset of H (0,1) 

such that 0 < a <. a(x) < b for each a € Q. (Here a 

and b are given bounds 7o r elements in Q.) 

5 N 

Turning to approximations for (6), we let S ( A ) 

denote the set of quintic splines (C 4 functions that 
are piecewise polynomials of degree 5 — see [6,10] 
for notation) corresponding to the partition 

i N = 1x.} N , X. ■ j/N. of [0,1]. We then define 

J-j-0 3 

Z (a) as the subspace of Z(a) given by 

Z N * Sg(i N ) x Sg(fl N ) where S*(a N ) s (s e S 5 (a N )|s(n) » 

0 2 s(n) * 0 , n * 0.1). Let P N (a) : Z - Z N be the 

u 

orthogonal projection of Z(a) onto Z (a) and define 
approximations to A by A N (a) = P N (a)A(a)P^(a). The 
approximations to (6) in Z (a) are then taken to be 


iV) = A N (a)z N (t) + P N F(t) t > 0, 

Z N (0) - P N Zg • 


(8) 


The approximate estimation problem Is that of mini- 
mi z i ng 

J N ( a ) - l |z'J(t i ;«)(x J ) - P^jl 2 (9) 

1 * j 

over Q. 

In order to establish a convergence theory for 
such Galerkln-type estimation schemes, it Is not dif- 
ficult (see [5,8,9]) to argue that It suffices to show 

that ^(t;a k ) — ► z(t;a*}, N ♦ k -► -, for any 
k * 

sequence {a } In Q converging to a In Q. The topology 
for this convergence depends on the problem at hand 2 
and in the present case. It is desirable to use the H* 
topology on Q. Using fundamental estimates from spline 
approximation theory (in particular, minor modifica- 
tions of estimates given In Lemmas 1.18 and 1.19 of 
[11] suffice) and fundamental properties of dissipative 
operators, one can readily employ the Trotter-Kato 
theorem (see [5]) to establish the desired convergence 


z N (t;a k ) -► z(t;o*) whenever a k -► a* In Q. Further- 
more, letting o N be a solution of minimizing (9) over 

Q, If we have (or a subsequence a t) converging to 
some a in Q (this requires a compactness property for 
Q), It can be argued that a Is a solution to the 
original parameter estimation probltm involving (6) and 

(7). Hence assuming that Q Is compact In H 2 , we have 
a sequence of problems which approximate In a desired 
sense our original problem. 


As we have already observed, the problem of 
minimizing (9) over Q subject to (8) Involves mini- 
mization over an Infinite dimensional parameter space 
and a further approximation Is necessary to obtain an 
Implementable computational procedure. For this we 
use cubic Interpol atory spline approximations to the 

elements of Q. For a given partition A* 




y, * j/M, of [0.1]. let I be the cubic spline inter- 
J M 

polatlon map (see [12, p.48]) corresponding to A and 

define Q M 5 I M (Q). Then for a given function a €Q, 


I n (a) involves the values a(yj)* j 
Da ( 0 ) , Da( 1 ) , SO that I M : Q 


C 1 topology on Q and the C° topology on Q M 


0,1,. ..,N, and 
Q M Is continuous In the 
If Q Is 

compact In the H* topology, then It Is therefore easily 
M 

seen that Q Is conpact in C(0,1) end furthermore has 
a representation 

M+3 


nM 


a: [0,1] ♦ R^l a 


nrg II 

l 

j-1 3 3 




where {Bj| are the cardinal cubic basis elements and 

the are compact subsets of r\ It follows that 0 M 
J 2 

Is compact In the H topology and hence the problem of 

u 

minimizing J over Q can be readily replaced by the 
finite dimensional state space - finite dimensional 

parameter space problem of minimizing J N over Q M . 

Under additional smoothness assumptions on the elements 
of Q, one can employ standard estimates for inter- 
polator splines to argue that Q M approximates Q in an 

appropriate sense and that solutions of the problem for 
N M 

J over Q approximate those of the problem of mini- 
mizing J over Q. More precisely one can establish the 
following convergence results. 

Theorem . Suppose Q C H 4 (0,1 ) with j|D 4 a| 0 } a< Q bounded 

and suppose further that 0 Is compact In the H 2 (0,1) 
topology. Let Q M - I M (Q) be the cubic spline Inter- 
polator approximations to Q and let, for each N and M, 
a N,M denote a solution of the finite state-finite 
parameter gpace problem of minimizing J N of (9) over q* 

Let U k * j l be any convergent (In the H 2 topology) 
subsequence with limit a*; then a* Is a solution to the 
problem of minimizing 0 of (7) over Q. 

Numerical Results 


We have developed and tested software packages 
based on the qulntlc/cublc spline scheme discussed 
above. The packages were modifications of those 
(described in some detail in [6]) we have used for con- 
stant parameter estimation problems (e.g., the IMSL 


f 


2 



package OSCAR Is used to solve the ordinary differen- 
tial equations (8), while the IMSl version - ZXSSQ • 
of the Levenberg-Narquardt algorithm Is used In mini- 
mizing J N over </V To test the method and software, 
"data" was generated In the following manner: A 
"true” parameter function « was chosen and an 
Independent numerical algorithm (finite differences) 
was used to generate corresponding solution values 

y 1J * 4 subs#t of "Mch was then used as 

observations or data In our Inverse algorithms. Me 
present here the results from one of our test examples. 
Additional numerical results along with a complete 
theoretical treatment will be given In a paper cur- 
rently In preparation. 

Example . We take (2) with «*(x) • .15 + .10 
tanh(5(x-.5)), a constant load f ■ 10. and assume that 
the beam Is Initially at rest (,■,■0). Observations 
at points Xj * .1..2......9 and times tj ■ .1..2,..., 

1.0. were used In the least squares criterion (9). 
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